Abstract 



There is growing evidence in favour of the temporal-coding hypothe- 
sis that temporal correlation of neuronal discharges may serve to bind dis- 
tributed neuronal activity into unique representations and, in particular, 
that 8 (3.5-7.5 Hz) and 5 (0.5 <3.5 Hz) oscillations facilitate information 
coding. The 9 and 5 rhythms are shown to be involved in various sleep 
stages, and during anaesthesia, and they undergo changes with the depth of 
anaesthesia. We introduce a thalamocortical model of interacting neuronal 
ensembles to describe phase relationships between and 5 oscillations, espe- 
cially during deep and light anaesthesia. Asymmetric and long range interac- 
tions among the thalamocortical neuronal oscillators are taken into account. 
The model results are compared with the experimental observations of Mu- 
sizza et al. J. Physiol. (London) 2007 580:315-326. The 5 and 6 activities 
are found to be separately generated and are governed by the thalamus 
and cortex respectively. Changes in the degree of intra-ensemble and inter- 
ensemble synchrony imply that the neuronal ensembles inhibit information 
coding during deep anaesthesia and facilitate it during light anaesthesia. 
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1 Introduction 



Neuronal communication and synchronization are crucially important fea- 
tures of the co-operative interaction between neuronal ensembles, endowing 
the brain with the marvelous capability known as cognition. The time scales 
of human motor and cognitive events are comparable to that of the EEG 
field dynamics arising from synchronous activity, observable from within or 
outside the brain (Jlj). Since EEG signals correspond to the averaged activ- 
ity of large cell populations, their fluctuations can support cognitive activity 
only if they are sufficiently synchronized: effective communication of neu- 
ronal ensembles can be achieved only if they are oscillating in a synchronized 
manner. Such behavior tags those neurons to represent particular cognitive 
tasks, e.g. in relation to a perceptual object (0). 

In this article we present a model of interacting thalamocortical neuronal 
ensembles in an attempt to account for the behaviour of 5 and 9 waves during 
anaesthesia. Our model is motivated by an attempt to tackle the problem of 
anaesthetic awareness, to provide the basic understanding needed to prevent 
people inadvertently awakening during surgery. We therefore focus on ways 
to identify and characterize the states of deep and light anaesthesia. Our 
starting point model is the experimental observation that both the cardio- 
respiratory and respiratory-5 interactions change with depth of anaesthesia in 
rats ((3)0) following the administration of a single bolus of ketamine-xylazine. 
During the ensuing deep phase of anaesthesia (~45 min) the amplitude of 
5- waves is strongly pronounced. As the subsequent light phase of anaesthesia 
(~25 min) is entered, the 5-waves disappear and the amplitude of 0-waves 
increases. We use the phase dynamics approach to propose a model that 
incorporates regional ensembles of neurons that are connected both among 
and between themselves. 

Further, we investigate the role played by the neuronal ensembles in 
temporal coding during deep and light anaesthesia from the model results. 
Neuronal synchrony is often associated with an oscillatory pattern of sig- 
nals (oscillation-based synchrony) . The frequencies of such signals generally 
cover a broad range and, more importantly, exhibit a marked state depen- 
dence. In other words, synchrony and temporal coding of information in 
selective frequency bands are two different sides of the same coin. We use 
this idea and establish the link between the conditions under which tempo- 
ral coding is believed to occur (synchrony) and the level of arousal of the 
brain. 
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1.1 Physiological background 

During periods of slow wave sleep/anaesthesia, widespread synchronized os- 
cillations occur throughout the thalamus and the cortex. The slow oscilla- 
tions (< 1Hz) that occur during natural sleep and ketamine-xylazine anaes- 
thesia are an emergent network property of neocortical neurons involving 
the TC and RE neurons (J5|). The clock-like 5 (1 — 4Hz) are produced in 
the thalamus and are strongly synchronized by the RE and the thalamocor- 
tical volleys (@). A small proportion of the TC neurons also display slow 
and 5 oscillations (0). Thus the high amplitude, low frequency 5 and slow 
oscillations are found to be associated with highly coherent activities of the 
cortical, RE and TC neurons. The dynamics of 5 and 9 waves can therefore 
be represented by interacting ensembles of neuronal oscillators. 

The physiological circuitry is as shown in Fig. [1] (left) (0)- The pyra- 
midal (PY) neurons of the cerebral cortex are connected among themselves 
(excitatory) by the interneurons (IN) which also form components of the 
cortical ensemble. The TC neurons receive sensory inputs which are relayed 
to the appropriate area of the cortex through ascending thalamocortical fi- 
bres (indicated by the upward arrow). The RE neurons wrap most of the 
dorsal and ventral aspects of thalamus (j3) and act as bridge between the 
TC neurons and the thalamic neurons. The dense axons of the RE neurons 
innervate the TC neurons (0). The corticothalamic fibers (indicated by the 
downward arrow) also leave collaterals within the RE nucleus and dorsal 
thalamus. The RE neurons thus form a network that surrounds the tha- 
lamus. It receives a copy of nearly all thalamocortical and corticothalamic 
activity, and projects connections solely to neurons in the TC region (0). In 
turn, the axons of the TC neurons give rise to collaterals in the RE nucleus 
while the parent axon passes through the cerebral cortex (|10l). 

2 The model 

Based on the physiological phenomena that take place during anaesthesia, 
as well as on anatomy, the model system considered is shown schematically 
in Fig. [1] (right). The oscillators in the three ensembles, namely the cortical, 
the thalamic reticular and the thalamocortical relay neurons, have different 
mean natural frequencies and their interactions are characterized by intra- 
population and inter-population coupling parameters. 

Each neuron in the ensembles is considered as an oscillator whose mem- 
brane potential is the oscillating variable. The couplings represent the 
synaptic connections between them. We reduce the system to a phase model, 
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one of the simplest yet accurate models for weakly coupled nonlinear oscil- 
lators, with the coupling being introduced through the phases. In doing so, 
we make use of the fact that there exists a degree of coherence between the 
membrane potential oscillations and the action potential firings. Although 
they are not equivalent, the action potentials are triggered at a certain phase 
of the membrane potential oscillations. This reasoning leads to the following 
set of equations 
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where 9 ( . 1 ' 2 ' 3 ^ are the phases of the ith oscillator in the C, TC and RE 
ensembles, respectively, and N refers to the ensemble sizes. The parameters 
A c , A tc and A re quantify intra-ensemble couplings and B c , B tc , Cf C and B re 
quantify the inter-ensemble couplings. The natural oscillator frequencies 
u;^ 1 ' 2 ' 3 ^ are assumed to be distributed with central frequencies Co 1 ' 2,3 . The 
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(77) ' ' ' ' (£)) = 2Z)( li 2 i 3)5(t — ^(1,2,3) are the noise intensities. 

For convenience one can define complex- valued, mean field, order parameters 

as r (i,2,3) e ^ (1 ' 2 ' 3) = jj Ylf=i e * j ■ Here V'1,2,3 (t) are the average phases of 
the oscillators in the respective ensembles and 7^2,3 (t) are measures of the 
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coherence of the oscillator ensembles, which vary from to 1. With these 
definitions, Eqs. (H|) become 





= 4 1 * 


— rij4 c sin(( 




- V'l + a) 






—r<iB c sin(6 


iw . 


- ^2 + a) + ?7i , 


(2) 

i 


= k (2 » 


- r 2 vli C sin( 


if 


- -02 + a) 






-T\B tc sin( 


of 


- + a) 






-r 3 C ic sin( 


f» 


/2 s ) 

- ^ 3 + a) + ?7> 


(3) 
i 




— r%A re sin 


,,(3) 


- 1P3 + a) 






—r2B re sin( 


Of 


- ^2 + a) + »7i 



(2) 



2.1 Numerical methods 

A fourth-order Runge-Kutta routine is used for the numerical simulation 
with the initial phases equally distributed within [0, 2ir]. The results are 
normalized to "real time" and the system is simulated for the equivalent of 
60 mins with N = 10000. We will refer to synchrony within an ensemble as 
intra-ensemble synchrony, and that between ensembles as inter-ensemble 
synchrony. The amount of intra-ensemble synchrony is measured by the 
mean field parameters: ^(1,2,3) = implies that there is no synchrony in 
the corresponding ensemble; ^(1,2, 3) = 1 indicates complete synchrony; and 
< ^(1,2,3) < 1 corresponds to partial synchrony. The greater the value 



of ^(1,2,3)1 the more oscillators are oscillating in synchrony (I 111 . Il2l ). On 
the other hand, inter-ensemble synchrony occurs when oscillators from two 
different ensembles entrain to a frequency/phase window. In general, inter- 
ensemble synchrony can be quantified by a constant difference in the mean 
phases ^(1,2,3)- However, this measurement will be valid only if both the 
ensembles are completely locked to each other. For partial synchrony, the 
difference in the mean phase will be oscillating, showing no synchrony. Hence 
we identify inter-ensemble synchrony by looking at the time evolution of the 
ensemble-averaged frequencies. Also, for increasing inter-ensemble coupling 
par ameters, a decrease in rn 23) is a signature of inter-ensemble synchrony 
0). 

The inter-ensemble and intra-ensemble coupling parameters are swept 
linearly in time in order to mimic the effect of decreasing concentration of 
anassthetic agent. This corresponds to the ability of anassthetics to affect 
thalamocortical signalling, as is well recognized from in vivo electrophysio- 
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logical work on animals ()14l ). Thus, during deep anaesthesia, the neuronal 
oscillators are very weakly coupled so that they do not interact much (due 
to the anaesthetics blocking the signalling pathways). As the anaesthetic 
concentration wears off the signalling pathways become unblocked and the 
couplings between the neuronal oscillators and ensembles become stronger. 
Hence we can simulate the transition from the deeply to the lightly anaes- 
thetized state by starting with a very weak coupling strength (characterizing 
deep anaesthesia) and continuously sweeping them to stronger values with 
time. The starting values of the coupling parameters are given in Fig. [2J 
We introduce asymmetry in the model by a phase shift < a < yr/2 . 
For the values of the mean frequencies we follow Amzica et al. (la, lla ): 
= 3 Hz, ci>( 2 ) = 1.5 Hz and a)^ = 1 Hz, with Lorentzian distributions 

Our choice of values for the mean frequencies reflects the intrinsic prop- 
erties of the neurons. The cortical neurons are endowed with intrinsic prop- 
erties which could be reflected in field potential recordings as S activities. 
The frequencies of these faster oscillations evolve around the upper limit of 



the delta band (mainly 3-4 Hz) (|15l . Ilq . Il7l ). The clock-like 5 oscillations 



(1-4 Hz) are mainly produced in the RE neurons ((B). Individual TC neu- 
rons are capable of producing Ca 2+ spikes and associated action potentials 
at frequencies of 0.5-4 Hz (|6|). 



3 Results 

The time evolution of the mean frequency of each ensemble is plotted in Fig. 
[2] right. For comparison, the experimental 5 and 9 frequencies, obtained by 
wavelet analysis of the EEG signals from anaesthetised rats (|3|), are plotted 
in Fig. [2] (left). In the latter experiments, the 5 activity was found initially 
to be of higher amplitude compared to that of the 9 waves. In particular, 
during the deep phase of anaesthesia, there occurs a strong, high amplitude, 
5 activity which greatly diminishes on entry to the light phase (at ~ 45 
min). The 9 activity runs independently throughout the whole period of 
anaesthesia, but at a much lower amplitude compared to that of the 5 waves. 
Note the differences in the amplitudes corresponding to the color codes in 
Fig. [2] left, which makes the 6 activity visible in the top panel. 

The model results suggest that it is the TC and RE ensembles that gen- 
erate the 8 activity during deep anaesthesia. This 5 is of high amplitude due 
to the strong synchrony in and between the TC and RE ensembles. Syn- 
chrony between the TC and RE ensembles can be easily established even 



5 



with a smaller value of inter-ensemble coupling, because the frequency dif- 
ference between them is relatively small. Fig. [3] shows the time evolutions 
of the mean field parameters ^(1,2,3) °f t ne three ensembles corresponding 
to the simulated frequencies. During deep anaesthesia, the fractions of os- 
cillators that oscillate in synchrony in the TC and RE ensembles are higher 
than in the C ensemble. This results in a strongly pronounced 5 while the 
amplitude of 9 remains low due to very low number of oscillators oscillating 
in synchrony in the C ensemble (see Fig. [3|). 

On the other hand, for the thalamus (TC+RE) to be synchronized with 
the cortex, a relatively higher value of coupling strength is required. This 
occurs at ~ 45 min, at which point the TC and RE frequencies shift to 
a higher value to join the C ensemble and produce 9 activity. This is the 
reason for the sudden diminution of the 5 and appearance of 9 at ~ 45 min. 
Of course, the exact time of occurrence depends on the value of the coupling 
parameter that we choose. As a consequence of the shift (inter-ensemble 
synchrony), the amount of intra-ensemble synchrony is reduced in the TC 
and RE ensembles (fl3l ) and hence the 9 appears as a low-amplitude activity. 
We thus quantify the amplitudes of the characteristic 5 and 9 activities by 
the mean field parameters ^"(1,2, 3) which reveals the amount of synchrony 
in each ensemble and hence the amplitude of the oscillations. The intra- 
ensemble and inter-ensemble synchronisation mechanisms discussed here are 
generic to systems of coupled oscillator ensembles (U, H, 13, H, 19, 20, 21 



23,123)- 



3.1 Phase coding and depth of anaesthesia 

The 5 and 9 phases play significant roles in coding information when the 
brain is not in an aroused state, characterized by desynchronized EEC The 
level with which the phases are synchronized determines the ability to code 
information. However, it is the inter-ensemble synchrony that plays the 
crucial role in temporal binding (rather than the intra-ensemble synchrony). 
Since synchrony is supposed to enhance the significance of responses of the 
neurons (j2), it is obvious that synchronized discharges will have a stronger 
impact than temporally disorganized ones (24). Thus the 9 phase is found to 
play a crucial role in inhibiting information by being poorly synchronized in 
the cortex. More importantly though, the cortex and the thalamus are not in 
synchrony with each other (see Fig. [3]). Due to this, binding of information 
cannot be achieved, which is why consciousness and cognition are absent 
during the deep phase of anaesthesia. The strong 5 waves keep the cortex 
and the thalamus out of phase with each other during deep anaesthesia. 
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On the other hand, the thalamus enters into synchrony with the cortex on 
emergence from anaesthesia. This corresponds to the state of awareness 
when information can be successfully coded due to the emergence of inter- 
ensemble synchrony, characterized by the 9 phase (see Figs. [2] and [3]). 



4 Discussion 

The model results indicate that the 5 and the 9 activities observed experi- 
mentally are separate in terms of their generation and frequency. The 5 and 
9 activities are found to occur mainly in 0.5 — 3.5 Hz and 3.5 — 7.5 Hz bands 
in agreement with recent observations (0) and earlier reports ( 15. EU). The 



dramatic diminution of the 5 amplitude and the simultaneous appearance 
of 9 activity characterize the transition from deep to light anaesthesia. 

In the experiment (0), the amplitudes of the 5 and 9 waves differed by 
more than a factor of 10. In order to reveal 9 (which otherwise would have 
been lost in the noise level of 5) two separate figures were therefore plotted 
with different amplitude scales. Our model results suggest that 9 activity 
is present during both deep and light anaesthesia. However, during deep 
anaesthesia the 5 activity is highly synchronized, whereas 9 activity is poorly 
synchronized (see Fig. [3j). The effect on their relative amplitudes is that 
the 9 activity cannot be seen during the deep anaesthesia in the experiment, 
consistent with the observations, despite its presence as revealed by the 
model. 

Although varying the coupling parameters to mimic the effect of de- 
creasing concentration of the anaesthetic agent helps us to understand and 
compare the model results with reality, it introduces oscillations in the fre- 
quencies. As is evident from Fig. [21 the mean frequencies of each of the 
ensembles undergo noisy oscillations, while remaining within the 5/9 bands. 
The effect of varying coupling parameters affects the frequencies directly be- 
cause the model equations are for the oscillator phases and the frequencies 
are calculated from the phases themselves. Likewise, the oscillations seen 
in the order parameters (Fig. [3]) are introduced mainly by the change in 
inter-ensemble coupling parameters. 

Since the model considers only the phases, and not the amplitudes, we 
are unable to make a direct comparison with the experimental results. How- 
ever, we make use of the theory of synchronization to quantify the amplitude 
in terms of the amount of synchrony in each of the ensembles. That is, the 
more the synchrony, the more oscillators are oscillating in phase with each 
other, and hence the higher the corresponding oscillation amplitude. On the 
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other hand, the model does offer the advantage that we are able to identify 
those neuronal groups that are responsible for the generation of 5 and 9 
activity during anaesthesia. 

Although we do not challenge the importance of the microscopic details 
underlying general anaesthesia and the Hodgkin-Huxley formalisms (|25l ). the 
model results support the hypothesis that consideration of macroscopic dy- 
namics with asymmetry is important. We therefore suggest that models of 
similar kinds can be used to explain experimental observations, not only in 
anaesthesia, as here, but also in other cognitive and behavioral states that 
involve huge numbers of neurons functioning in groups. Detailed work has 
been done in the field of neuronal mass modeling and a wide range of EEG 
generative models have been proposed 



30, I3J, Our 
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model is capable of generating the wide range of EEG oscillatory behav- 
iors reported earlier (33, 34, 35, 3^, 37). Our findings uniquely identify the 
link between neuronal synchrony and temporal coding especially in terms of 
inter-ensemble synchrony. An additional advantage of the model is that it 
yields insight into the synchronization mechanisms underlying the 5 and 9 



waves. 



5 Conclusion 

In summary, we have introduced a model involving asymmetrically inter- 
acting ensembles of C, TC and RE neurons in order to understand the 
mechanisms underlying the generation of 5 and 9 waves during anaesthesia. 
The model results are compared with those from experiments (J3j). The TC 
and the RE ensembles were found to be responsible for the generation of 
high amplitude 5 waves during deep anaesthesia. The C ensemble is engaged 
with the 9 activity. The transition from deep to light anaesthesia is found to 
be marked by a frequency shift in the TC and RE ensembles, caused by the 
increase in the coupling strengths. Also, the 9 activity is found not to be 
as strongly synchronized as the 5 activity. The similarities and differences 
between the model results and the experimental results were discussed. Fur- 
thermore, the model illuminates the phenomenon of temporal coding of in- 
formation, in particular by the 5 and 9 frequencies. It reveals the role played 
by the neuronal phases in inhibiting information during deep anaesthesia and 
coding the sensory information during light anaesthesia. Although our main 
motivation for introducing the model was to understand the mechanisms 
giving rise to the generation of 5 and 9 waves during ketamine-xylazine 
anaesthesia in rats H, 0) , it can also be used to elucidate the mechanisms in- 
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volved in the generation of brain waves for other anaesthetics and also during 
slow wave sleep. The results derived from the model will have implications 
for the understanding of fMRI and MEG dynamics. 

The study was supported by the EC FP6 NEST-Pathfinder project 
BRACCIA and in part by the Slovenian Research Agency. 
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Figure Legends 



Figure [TJ, 

Left - The arrangements and connectivity of four types of cells, namely 
the TC, RE, C (comprised of the IN and PY cells). Pre refers to external 
pre-thalamic afferent sensory inputs. The upward and downward arrows 
represent the thalamocortical and corticothalamic fibres respectively. Af- 
ter (0), with permission. Right - Schematic representation of the model. 
There are three ensembles: cortical neurons(C), thalamocortical relay neu- 
rons (TC), and thalamic reticular neurons (RE). They interact both among 
themselves (represented by the circle in each case) and between each other 
(represented by the arrows). The upward and the downward arrows indi- 
cate the thalamocortical and corticothalamic connections. External afferent 
(pre-thalamic) inputs to the thalamus are denoted as Pre. 

Figure [2J 

Left - The time evolution of the characteristic EEG 5 and 8 frequencies 
during Ketamine-xylazine anaesthesia, analyzed by wavelet transform. Re- 
produced from (0), with permission. Right - The time evolution of the 
characteristic 5 and 8 frequencies displayed by three ensembles C (blue), TC 
(red) and RE (green), as obtained from the model. The (starting) values of 
the (coupling) parameters were A c = 0.8, B c = 1.2, A tc = 0.9, B tc = 0.45, 
Cu. = 0.9, A re = 0.2, B re = 0.65, a = 0.9, D 1 = 0.1, D 2 = 0.2, D 3 = 0.15 
and 7 = 0.4. 

Figure |3J, 

Mean field parameters rn 2,3) of the three ensembles, plotted as functions of 
time, corresponding to the frequencies plotted in Fig. [2] (right) for the same 
values of parameters. 
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Figure 1: 
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